The Arbitrary Trajectory Quantization Method 
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The arbitrary trajectory quantization method (ATQM) is a time dependent approach to quasi- 
classical quantization based on the approximate dual relationship that exists between the quantum 
energy spectra and classical periodic orbits. It has recently been shown however, that, for polygonal 
billiards, the periodicity criterion must be relaxed to include closed almost-periodic (CAP) orbit 
families in this relationship. In light of this result, we reinvestigate the ATQM and show that at 
finite energies, a smoothened quasiclassical kernel corresponds to the modified formula that includes 
CAP families while the S function kernel corresponding to the periodic orbit formula is recovered as 
E — > oo. Several clarifications are also provided. 
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I. INTRODUCTION 



Semiclassical quantization methods developed over the 
past three decades generally rely on the energy domain 
approach. The spectral density, p(E) = J2 n $(■& ~ E n ) is 
thus expressed in terms of periodic orbits of the underly- 
ing classical dynamics and this is commonly referred to 
as the Gutzwiller trace formula ["J. The crossover from 
the time domain to the energy domain and the trace for- 
mula necessitates several approximations in the form of 
stationary phase integrations and one of these selects pe- 
riodic orbits as the sole classical ingredient. The recipe 
works well in most cases since closed non-periodic orbits 
contribute with a lesser weight and can be included only 
as a correction. There are instances however when closed 
non-periodic orbits contribute with weights comparable 
to periodic orbits and the "periodicity criterion" must 
then be relaxed to accommodate them [0,0. Polygonal 
billiards provide such an example and perhaps hold the 
key to generic intermittent behaviour. In these systems, 
a slight change in the internal angles results in the de- 
struction of periodic orbit families thereby giving rise to 
closed almost-periodic (CAP) families of orbits [0. Since 
this is true for each periodic family in every neighbour- 
ing polygon, CAP families actually outnumber periodic 
families of any given polygon. Further, they contribute 
with weights comparable to periodic families in the semi- 
classical trace formula and are hence indispensable for 
semiclassical quantization. This makes the energy do- 
main approach rather cumbersome to implement and a 
time domain approach is thus preferable. 

For polygonal billiards, there are two available ap- 
proaches in the time domain. The first of these devel- 
oped by Heller and Tomsovic J5| relies on a semiclassical 
construction of the time dependent propagator to eval- 
uate the auto-correlation function, C(t) —< ip\i/)(t) > 
where ^(q,t) =< q|^(t) >= £„ c n M^~ lEnt/h - 
f K s . c (q,q! ,t)ip(q! ,0)dq! . Here {</>„}, {E n } and {c n } are 
the eigenfunctions, eigenvalues and expansion coefficients 
respectively while _K" s . c (q, q',f) is the semiclassical prop- 
agator constructed using classical trajectories joining a 



pair of points (q, q') at each time t. A power spectrum of 
"0(q, t) thus yields the quantum eigenvalues. The method 
has been successfully applied to the stadium billiard but 
to the best of our knowledge, it has not been used for 
quantizing polygonal billiards. 

The second approach ||[7J is simpler to adopt and re- 
lies on classical propagation (or its quasiclassical Q adap- 
tation where necessary). In case of a polygon, it involves 
shooting arbitrary trajectories in various directions from 
a point (call it q' - see fig. [y, left) and at each time 
step, recording the (weighted) fraction, F{t) of trajecto- 
ries that are in an e neighbourhood of a point q. The 
peaks in the power spectrum of F(t) are then related to 
the quantum eigenvalues. 





FIG. 1. The arbitrary trajectory quantization recipe. In 
case of the polygon, trajectories with different linear momen- 
tum span the entire constant energy surface. In case of the 
circle, this is achieved by shooting trajectories with different 
angular momentum. In both cases, the (weighted) fraction of 
trajectories in the cell around q is recorded. 

Apart from the simplicity, the arbitrary trajectory 
quantization method (ATQM) (see ||) is perhaps the 
only semiclassical scheme that has been successful in de- 
termining the first few quantal energy levels of generic 
polygonal billiards M. However, the theory as presented 
in |6J leaves a number of questions unanswered. For one, 
the role of the e neighbourhood is unclear. Besides, as 
the theory uses periodic orbit quantization to relate the 
peaks in the power spectrum of F(t) with the quantum 
eigenenergies, the role of CAP trajectories must be clar- 
ified. We shall thus reinvestigate ATQM from this view- 



point and understand why it works. 

The plan of the paper is as follows. In section ||, 
we shall recapitulate the existing work on the arbitrary 
trajectory quantization method. The modifications that 
we shall carry out to account for the inclusion of closed 
almost-periodic orbits can be found in section III and 
this constitutes the main part of this paper. Finally, a 
discussion on errors and a summary of our results can be 



found in section IV 



II. ARBITRARY TRAJECTORY 
QUANTIZATION : A REVIEW 

The arbitrary trajectory quantization method relies 
on a suitable (quasiclassical) adaptation of the classi- 
cal evolution operator that propagates a density under 
a flow that we denote by q l . For polygonal billiards, the 
flow occurs on an invariant surface that is 2-dimensional 
and characterized by the two constants of motion, {E, tp} 
where E is the energy and (p denotes the second constant. 
The surface has the topology of a sphere with g (called 
the genus) holes where g can be determined from the in- 
ternal angles of the polygon. The motion on the invariant 
surface can alternately be viewed on a singly connected 
surface obtained by executing 2g cuts and with edges 
appropriately identified. 
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FIG. 2. The singly connected region for an L-shaped bil- 
liard consists of four copies with edges appropriately identi- 
fied. A trajectory originating near the 37r/2 vertex in 1 is 
plotted in configuration space using bold lines and the corre- 
sponding unfolded trajectory is also shown. The latter con- 
sists of parallel segments and the trajectory can be parame- 
terized by the angle ip that it makes for example with the q\ 
axis. 

As a trivial example, consider the rectangular billiard. 
The singly connected region is a larger rectangle consist- 
ing of four copies corresponding to the four directions 
that a trajectory can have and these can be glued ap- 
propriately to form a torus. As a non-trivial example, 
consider the L-shaped billiard of fig. g which is pseudo- 



integrable with its invariant surface having, g = 2. Al- 
ternately, the surface can be represented by a singly con- 
nected region in the plane and consists of four copies 
corresponding to the four possible directions an orbit can 
have and these are glued appropriately. A trajectory in 
phase space thus consists of parallel segments at an an- 
gle ip [QOl measured for example with respect to one of 
the sides. It will be useful to note at this point that the 
same trajectory can also be represented by parallel seg- 
ments at angles n — (p, n + ip and 27r — p. In general, the 
number of directions for representing a trajectory equals 
the number of copies, N, that constitute the invariant 
surface. 

Before considering the question of quasiclassical quan- 
tization, we first introduce the appropriate classical evo- 
lution operator. For integrable systems, this is easily 
defined as 



£*o<; 



2) 



C$1 C$2 6(0i - 0f)<5(0 2 
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(1) 



where 9\ and 62 are the angular coordinates on the torus 
and evolve in time as 6\ — u>i(Ii,l2)t + 6i with Wj = 
dH(I u I 2 )/dIi and h = ^£,.p.dq. Here T 4 ,i = 1,2 
refer to the two irreducible circuits on the torus and p is 
the momentum conjugate to the coordinate q. 

On a general two dimensional invariant surface param- 
eterized by ip, the classical propagator is expressed as H 



£*(^)ocA(q) = J d4 %- q*(¥>)) 0(q') 



(2) 



where q refers to the position in the singly connected 
region (see fig. ||) and q n (ip) is the time evolution pa- 
rameterized by ip as described above. We denote by 
{A n (t\ip)}, the eigenvalues of £'((/?) and from its mul- 
tiplicative nature, it follows that A n (i; p) — e x "^*'. We 
are interested here in the form of X n (<fi) and to this end 
we shall evaluate 



/ dip Tr C\<p) = / dip^2< 



,A„( v )t 



dip / dq <5(q - q*(<^)). 



(3) 



The delta function kernel ensures that the only orbits 
that contribute are ones that are periodic or close on 
the invariant surface after a time t. The q integra- 
tions arc thus simpler to perform if we transform to 
a local coordinate system with one component parallel 
to the trajectory and the other perpendicular. Thus 
S\\ (<?|| - qp = J2 P J2 r h S (t ~ tT p) wnere v is the velocity, 
T p is the period of the orbit and r is the repetition num- 
ber. Similarly, for an orbit of period rT p parameterized 
by the angle <p p , S ± (q ± -q r ± ") = 6(<p-(p p )/\dq ± /d(p\ v=ipp 
where \dqj_/dip\ — rl p for marginally unstable bil- 

liards. Putting these results together and noting that 
each periodic orbit occurs in general at N p different val- 
ues of tp, we finally have 



/^£ eA " wt = £E^(*-^) 



^E?^- r « ( 4 ) 
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H(I U I 2 ) = * 2 {Ii/LJ + Il/L 2 2 ) where L X ,L 2 are the 
lengths of the two sides. With I\ = \f~EL\ cos((p)/ir and 
I 2 = \Z~EL 2 sin(<£)/7r, it is easy to see that at a given en- 
ergy, E, each torus is parameterized by a particular value 
of <p. Thus 



where I = tv and the summation over p refers to all prim- 
itive periodic orbit families with length l p and occupying 
an area a p . Note that in Eq. (J|), we have replaced N p by 
N since for most long orbits, N p ~ N. Also, we have ne- 
glected the influence of isolated orbits to simply matters. 
We shall continue to make this approximation through 
the rest of this paper and justify its use at the end. 

In some cases, it is possible to interpret the periodic 
orbit sum in Eq. (0) starting with the semiclassical (E — > 
oo) trace formula for marginally stable systems : 



A ft" U)) — g i27r *v / ^(™i cos(ip)/L 1 +n 2 sin(tp)/L 2 ) 



(8) 



and 
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rfip). (5) 



Here p a v(E) refers to the average density of quantal 
eigenstates, k = \/E, l p is the length of a primitive pe- 
riodic orbit family, fi p = irn p + i/ p ir/2, n p the number of 
(phase-altering) bounces that it suffers at the boundary 
and v p the number of caustics encountered by the orbit. 
Note that in the Neumann case, n p — since there is 
no phase loss on reflection while for polygonal billiards, 
v p = 0. For convenience, we have U = 1, v = 1 and the 
mass to = 1/2. Starting with the function 



where I — 2ty/E, tan(/i n ) = n.i-L 2 /(n 2 £i) and E n = 
it 2 (n\/L\ + n^/Ll). On separating out n = (0,0) from 
the rest, restricting the summation to the first quadrant 
of the integer lattice and noting that for a rectangle bo = 
1/4, Eq. (|) follows. 

The classical evolution operator thus serves to deter- 
mine the Neumann spectrum in polygonal billiards. Ap- 
propriate modifications however need to be made for the 
Dirichlet spectrum or for systems that have caustics (the 
circle billiard is an example) and the construction of the 
evolution operator is then guided by the nature of the 
semiclassical trace formula (Eq. g) . 

The quasiclassical evolution operator, C qc , linking 
the classical eigenvalues to the the desired semiclassical 
eigenvalues can be defined as 



_ r°° £* (<p)od)(q)= [ da' 6(0-0'* (<p))e 

J A 



■in{t) , K— iv{t) -. 



(10) 



(6) 



where f(x) = J ■£- cos(x — tt/4) and < A < Eq, it is 
possible to show using Eq. (g) that for polygonal billiards 

ten 



J ^E* 



,A„(¥>)t 



27riVbo + 27rJV^/(v^rO- (?) 



where {E n } are the Neumann eigenvalues of the sys- 
tem and 60 is a constant 11 I3|. Thus \ n (p) — 
WEn sin((p) [O. This is the central result of M when 
the Maslov phases are zero. 

For integrable polygons, Eq. (Jh can in fact be de- 
rived directly starting from Eq. (h). The eigenfunctions, 



{4>n(^i>^2)}i on the torus are such that 



?i,< 



D 



A„(i)</>„(0i, ^2) where A n (t) = e l<Xn . On demanding that 
4> n (di,02) be a single valued function of (61, 62), it fol- 
lows that </>„(0i,0 2 ) = e^ ni9l+n20 ^ where n = (m,n 2 ) 
is a point on the integer lattice. Thus the eigenvalue, 
A n (t) = exp{££(niWi + n 2 uJ2)}- 

To illustrate the relationship between J dp^e^^ 1 
and {E n }, consider a rectangular billiard for which the 
Hamiltonian expressed in terms of the actions, Ji, J 2 is 



where v(t) = z/(q' (tp)) and n(t) — n(q[ (p)) count re- 
spectively the number of caustics and (phase altering) 
reflections encountered by the trajectory q' (ip) in time 
t. As before, the multiplicative nature of £„ c (ip) implies 
that its spectrum is of the form {e^^)'} and it remains 
to be shown that for the quasiclassical operator defined 
in Eq. (fro), {A„} has a one-to-one correspondence with 
the appropriate quantum spectrum. 
As before, we shall evaluate 

= f dip f dq6(q_-q t (ip)) e- in{t)n ' iu(t) ^ (11) 

The trajectories that contribute are once more periodic 
due to the delta function in the kernel. The q and <p 
integrations can be performed along similar lines and we 
finally have 



/oo 
^^-"EE^ 
p r=l P 



6(1 — rl p )e 



-irn v Tr—trv„ 



(12) 



Starting with the function J2 n g(y/E^l)exp(—(3E n ), it 
follows from Eq. (j|) that for (3 — > + , 

J dip Tr C\ c =Y,K(t) = 2^VC + 2vriV^ 9 (v / ^0 

(13) 

where {-En} now refers to the desired quantum spectrum, 
g(x) — y/2/(irx) exp(— ix + 7r/4) and C is a constant. 
Since 



9(x) 



-ix sin(tp) 



(!^p 



for large x, it follows that for v = 1, 
X n ((p) = zv / ^"sin(( / 3) 
Eq. (p"5|) forms the central result of [E 



(14) 
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FIG. 3. The power spectrum, S(k) for the triangular en- 
closure with base length unity. The top set of arrows mark 
the first ten Dirichlet eigenvalues of the triangle. The bottom 
set are obtained from Bogomolny's transfer operator method 
[|L4| . Note that the non-zero cell size (e) corresponds to the 
hat-function kernel rather than the delta function kernel used 
so far. For more details, see text. 

To demonstrate this result, we shall consider a trian- 
gle with base angles (ir/5, 37r/10) and evaluate the power 
spectrum, S(k), of the (phase weighted) fraction of tra- 
jectories averaged over 150 cells of radius e = 0.03. The 
(p integration is performed using 300 trajectories each of 
which has a length 2048 x 0.025. The result is plotted 
in fig. |3| The peak positions of S(k) approximate the 
exact quantum eigenvalues well. Note however that the 
cells have a non-zero value of e and the power spectrum 
depends sensitively on the choice of this quantity. In 
the following section therefore, we shall develop an ap- 
propriate theory that takes into account the non-zero e 
requirement. 



III. ARBITRARY TRAJECTORY 
QUANTIZATION : MODIFICATIONS 

As pointed out, the existing theory does not require 
the e neighbourhood depicted in fig. [y. All one needs is 
the expansion of the kernel 



ff(q,q») 



~0™(q)</> r *(q') e 



i\/E n sin (ip)t 



(16) 



where 4> n are the eigenfunctions of C qc . The power spec- 
trum, S(k), of J dip K(q, q n (ip)) then has peaks at \/E n . 
In practice, this implies that e must tend to zero since 
the kernel has a delta function. The results however get 
worse as e is reduced to zero. A non-zero cell size on the 
other hand gives satisfactory results but is equivalent to 
a smoothened kernel p5[ . There is thus a gap in the 
present understanding and it is necessary to take a fresh 
look at the theoretical framework. 

The trouble in effect lies with the semiclassical trace 
formula (Eq. (||)) which neglects the contributions of 
closed almost-periodic trajectories. In case of a poly- 
gon, these occur in families with the same symbolic dy- 
namics and across which the action varies slowly. To 
see this, consider an arbitrary triangle T. In its imme- 
diate neighbourhood (obtained by changing the angles 
slightly), there exists an infinity of triangles {T^}, each 
with a distinct periodic orbit spectrum but having the 
same symbol sequence (in which the sides are visited) for 
times that depend on the differences in angles. Assume 
now that there exists a periodic orbit corresponding to 
a sequence Sk for the triangle T^. Then, for all other 
triangles in its neighbourhood, this sequence contributes 
to the semiclassical trace formula, an amount (nearly) 
equal to the periodic orbit contribution of T^> provided 
TTWitpe << A |j]. Here Adi is the angle between the ini- 
tial and final momentum of the orbit, u>i is the transverse 
extent of the family and A is the de Broglie wavelength. 
Thus, corresponding to every periodic family in each of 
the triangles {T^}, there exists an almost-periodic fam- 
ily in the triangle T whose contribution is comparable 
to that of periodic orbit families in these neighbouring 
triangles. 

The semiclassical trace formula should thus incorpo- 
rate CAP orbits and the modified expression is || : 



p{E) ~ Pav (E) +J2 



x cos(kli — 7r/4) 



sm(kAdiWi/2) 
kABiWi/2 



(17) 



where the sum over i runs over closed almost-periodic 
and periodic orbit families of transverse extent uii , U is 
the average length of such a family (taken as the length 
of the orbit at the centre of the band) while Adi is the 
angle between the initial and final momentum of the or- 
bit. Note that at any finite k, there exists a CAP family 
for which the de Broglie wavelength, A >> TTWiAOi so 



that the family contributes to the modified trace formula 
with a weight comparable to that of periodic families 

(0(1/^/2)). 

The modified trace formula (Eq. [17]) however fails to 
relate the quantum eigenvalues with the eigenvalues of 
C qc since its 8 function kernel chooses only periodic or- 
bits. The kernel function must therefore be suitably 
smoothened to achieve such a correspondence. 

In order that there exists a direct relationship between 
the eigenvalues of a quasiclassical operator, C qc and the 
quantum eigenvalues, {E n }, the kernel function corre- 
sponding to Cq C must be [pi 



K qc (<i,q[,t) =K s (q,q!,t)e 



-in(t)7T — il^(t) -^ 



(18) 



where 



Ka(q,q',<) 



sin (k{q\\ - q'^ip)) sin [k(q± - q'xif)) 



tt(9|| - 3|i*(v)) AQ± - <7l'(^)) 



(19) 



Note that lim^^, K s (q, q', t) = S(q- q'*(y)). The sub- 
script s thus denotes the smoothened kernel. As before, 
the multiplicative nature of C (ip) implies that its spec- 
trum is of the form {e x ^ ip ' t } and we shall now show 
that for the quasiclassical operator defined in Eq. (18), 
{A^} has a one-to-one correspondence with the appropri- 
ate quantum spectrum. 

Note that in general, gy— gy*^) = li+q±A8i — l = Agy 
for the i th (CAP or periodic) family while q± — gj_*((y9) = 
liip — Agj_ so that 



/ dip dq K s (q, q, t) = / dipdq\\dq± 



sin(fcAgy) sm(kAq±) 
7rAgy irAq± 

(20) 



A few approximations are now in order to keep the deriva- 
tion simple. First, we shall replace Agy in the denomi- 
nator by its mean value (at q± = 0) so that Agy ~ I — li 
Thus 



I" 1 sin(fcAgy) sin(kAe i w i /2' J 

d Q± — 7X7 ~ 2 TlTn 6 ki l - l i) 



7rAgii 



kA9, 



(21) 



where 8^(1 — li) = s'm(k(l — li))/(ir(l — li)). We shall next 
consider k finite but sufficiently large so that 



/ 



dip 



sin(fcAgj_) 
7rAg_L 



(22) 



The <7|| integration needs no approximation and yields 
/ dq\\ = l\ where l\ refers to the length of the orbit (prim- 
itive when A9i = 0). Finally then, 



J dip TrC^- ^N^frSiil 



h)e 



— irii'K—tVA 



(23) 



where /3* = sin(k A8 iiu i / 2) /(kA9iWi/ 2) and a, = Wilf. 
For a polygon, m is even for both CAP and periodic 
families while v^ = 0. Thus as k — > oo, Eq. (E3h reduces 
to Eq. (0). As before, we have neglected the contribution 
of isolated orbits from both the classical and semiclassical 
trace formulae. 

The eigenvalues of £* e can be related to quantum 
eigenvalues using Eqns. ( p3| ) and (Qj) and we merely state 
the final result : 



K(<P) =«V^sin(v?) 



(24) 



This is no different from our earlier result but in practi- 
cal terms, the use of the smoothened kernel justifies the 
use of the e neighbourhood since K s contributes substan- 
tially only in a small neighbourhood around q, the size of 
which is determined by k. Note that the correspondence 
demands that the quantity k in K s be identified with 
\/Krnax where E max is the maximum energy that one is 
interested in. Thus k also determines the time increment 
in the evolution of the kernel. 

Undoubtedly, several other smoothened kernels are 
just as appropriate. For instance, if one evaluates the 
fraction of trajectories in an e neighbourhood of q, the 
hat function of width e and height 1/e is the appropriate 
kernel to use. The role of k is then replaced by -. On 
the other hand, if one chooses a Gaussian centred at q, 
the appropriate smoothened kernel is 



(g ll — g||*tvJ) 2 (ix-glV)) 2 



^ s G (q,q',£) 



V / 2 7 rcr 2 



V27R7 2 



(25) 



On performing the three (ip, qn, q±) integrals similarly, it 
is possible to show that 



d ^ K ? = HrV-\ 



1 f wjMi 

/; I 3 I 2 3 / 2 CT 



(t-hlz 



V27R7 2 

(26) 



which approximates Eq. (E3J) for orbits with A9i small 
provided 1/tr is identified with k. Thus any smoothened 
delta function ought to work reasonably well. 

Finally, a few remarks about the neglect of the iso- 
lated periodic orbits are in order. First, these are fi- 
nite in number and have a smaller contribution in both 
the (quasi)classical as well as the semiclassical trace for- 
mula compared to periodic and almost-periodic families. 
Then again, diffractive contributions |l7f can be of the 
same order as isolated periodic orbit contributions in the 
semiclassical trace formula. Thus, in order to be able to 
compare the semiclassical eigenvalues and the eigenvalues 
of Cqc, we have consistently neglected isolated orbits. 

We provide some numerical results now. We first 
evaluate the function G(t) = J dip K qc for the 
(1997r/1011,3l7r/103) triangle of unit base length us- 
ing about 10000 trajectories for the ip integration. To 
achieve smoothening, we use the Hanning window func- 
tion together with a Gaussian damping and find the 



power spectrum, S(k), of G(i)e~' 3t . This is the inten- 
sity weighted spectrum, J2 n \cj) n (q)\ 2 5 s (k — k n ) (where 
6 S is a smoothened delta function), for a single pair of 
points (q, q' = q). Thus, the number of distinct and 
unambiguous peaks in S(k) depends on the value of k 
and |</>„(q)| 2 so that by changing q, a different set of 
peaks may be generated (see fig. 0). In fig. Ha, the nine 
highest peak locations (those above the dashed line) are 
compared with the exact quantum eigenenergies while 
in fig. 0p, the highest five are compared. In all cases, 
the agreement is good and the slight difference between 
peak locations and the exact quantum levels is due to the 
semiclassical nature of the calculation. 




FIG. 4. The power spectrum S(k) of G(t)e _,3t as a func- 
tion of k is shown for two different sets of points (q, q' = q) 
in (a) and (b) respectively. The arrows in (a) mark the ex- 
act quantum eigenenergies corresponding to the nine tallest 
peaks while those in (b) mark the five tallest peaks. The 
dashed lines are used to mark the demarcate the tallest peaks 
in both cases. Here k = 100, (3 = 0.01 and the length of each 
trajectory is 20487r/fc. There exist 37 eigenvalues for k < 50. 

The intensity weighted spectrum of fig. |4| thus demon- 
strates that arbitrary trajectories can be used to extract 
information about the quantum eigenstates. Can this 
however be used as an effective quantization method to 
determine all eigenvalues ? Fig. [3] indicates that this may 



indeed be possible. A natural way of achieving this is by 
integrating the intensity weighted spectrum over q JL8). 
In practice, only a few q values (typically 10-20) provide 
information of about 80 % of the eigenvalues in a given 
energy range. 



IV. DISCUSSION AND CONCLUSIONS 

It is worth recalling the sources of error in the quanti- 
zation recipe presented here. The theoretical basis pre- 
sented in sections || and III clearly indicates that the 
eigenvalues obtained using this method are at best "semi- 
classical" in nature since the modified trace formula, 
which connects the eigenvalues of C qc with the quantum 
eigenvalues, is only approximate as higher order correc- 
tions (0(l/k) due to isolated and diffractive orbits) have 
been neglected. In cases when the corrections are zero 
(such as in the rectangular or equilateral billiards), the 
ATQM does give exact results. However, there are ex- 
amples of other integrable billiards (such as the circle 
billiard) when corrections to the semiclassical trace for- 
mula are non-zero and the ATQM gives only the EBK 
eigenvalues |1S| . 

There can however be a further source of error in the 
approximation N p ~ N (see the discussion after Eq. ||) 
especially in non-generic systems where a significant frac- 
tion of long periodic orbits access fewer momentum direc- 
tions than the permissible number. However, we believe 
that in generic situations, the approximation is fair. 

Finally, since a comparison of this method with peri- 
odic orbit theory is inevitable, it must be reiterated that 
the "usual" periodic orbit theory neglects closed almost- 
periodic orbits and hence cannot give correct results in 
generic situations. The modified periodic orbit theory 
(or trace formula) does include CAP orbits. However, 
due to difficulties in enumerating them, energy domain 
quantization is expected to be quite cumbersome while 
the ATQM scores well due to its simplicity. There are 
non-generic situations however where closed almost peri- 
odic trajectories do not contribute significantly and this 
occurs in systems where the number of momentum di- 
rections accessible is small and A6i is therefore always 
"large". In such situations, diffraction effects assume 
greater significance and the ATQM may not be as ef- 
fective. 

In summary, we have provided a simple algorithm for 
determining the intensity weighted semiclassical spec- 
trum of polygonal billiards and shown that a smoothened 
delta function kernel for the quasiclassical evolution of 
densities is necessary at finite energies to incorporate the 
effects of closed almost-periodic orbits. 
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